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Abstract 

We investigate the structural and magnetic origins of the unusual ultrafast second-harmonic- 
generation (SHG) response of femtosecond-laser-excited nickel oxide (NiO) previously attributed 
to oscillatory reorientation dynamics of the magnetic structure induced by d-d excitations. Using 
time-resolved x-ray diffraction from the (| | |) magnetic planes, we show that changes in the 
magnitude of the magnetic structure factor following ultrafast optical excitation are limited to 
A (Fm) / {Fm) = 1.5% in the first 30 ps. An extended investigation of the ultrafast SHG response 
reveals a strong dependence on wavelength as well as characteristic echoes, both of which give 
evidence for an acoustic origin of the dynamics. We therefore propose an alternative mechanism 
for the SHG response based on perturbations of the nonlinear susceptibility via optically induced 
strain in a spatially confined medium. In this model, the two observed oscillation periods can be 
understood as the times required for an acoustic strain wave to traverse one coherence length of 
the SHG process in either the collinear or anti-collinear geometries. 
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As materials for spintronic applications, antiferromagnets (AFMs) provide the advan¬ 
tageous property of potentially faster spin switching times compared to ferromagnet^®^. 
By exploiting the exchange bias interaction that arises on an interface to a ferromagnetic 
materiaP, AFMs could act as fast switches in future data storage devices. Due to its high 
Neel temperature = 523 K, the AFM NiO is a promising candidate for such an applica¬ 
tion provided that the spin order can indeed be quickly and efficiently switched. 

Controlling the antiferromagnetic order in NiO appeared to become feasible after it was 
reported that ultrafast optical excitation could lead to a change in the magnetic anisotropy 
potential, inducing 90° flips in a substantial fraction of the spin population and a rotation of 
the macroscopic order parameter alongside a sudden decrease in spin ordei®®. Time-resolved 
SHG studies on bulk NiO further suggested that the dynamics were highly sensitive to the 
intensity and duration of the excitation pulse, resulting in transient redirection of spins os¬ 
cillating in either of two distinct directions with frequencies of 1 and 55 GHz, respectiveljP. 
This model was supported by the agreement of the observed frequencies with the frequen¬ 
cies associated with differences in anisotropy energy to the easy axis ground state known 
from neutron diffractiorP. Similar mechanisms leading to spin reorientation after optical 
excitation have been described for AFM compounds exhibiting a net magnetizatioiP. These 
materials, however, show a strongly temperature-dependent magnetic anisotropy that is not 
seen in NiO. The observations of AFM dynamics in NiO have so far only been made us¬ 
ing resonantly enhanced SHG of i?Ff-type, involving a magnetic-dipole transition. This is 
an elegant but indirect measure of magnetism in NiO which poses signihcant experimental 
challenged. 

Similar optical excitation levels in absorbing solids also generate coherent strain waves 
with frequency components similar to some of those observed in the SHG measurement^. 
For NiO in particular, Bosco et ah observed oscillations in time-resolved linear reflectivity 
measurements and Takahara et ah recently remarked at the close resemblance between fre¬ 
quencies of these coherent oscillations and the higher frequency oscillations reported for the 
SHG respons d^^*^^ ! It remains, however, unclear how such coherent strain waves can ade¬ 
quately explain all features observed in SHG measurements. In particular the magnitude of 
the changes in SHG are signihcantly larger than those seen in reflectivity, and the depen¬ 
dence on fluence appears to be substantially different. Linear reflectivity measurements also 
show no evidence of the 1 GHz oscillations seen in SHG. 
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By performing additional measurements and simulations, here we demonstrate that the 
dynamics observed in the SHG response to ultrafast optical excitation can be explained 
self-consistently as the result of coherent strain propagation. For this we used time-resolved 
x-ray diffraction to directly measure the evolution of spin order after excitation in this system 
and performed extended SHG studies that validate and quantify the strain wave induced 
origin of the SHG dynamics in NiO. 


I. STRUCTURAL, MAGNETIC AND DOMAIN PROPERTIES IN NIO 

Above its Neel temperature, NiO is paramagnetic and is found in a rock-salt crystal 
structure. Below T/v, Ni^^ spins align ferromagnetically within (111) planes, pointing along 
one of three equivalent {112} directions. The magnetization direction of neighboring (111) 
planes is antiparallel. Magnetostriction leads to small structural contractions along (111), 
giving rise to four possible structural domain orientations derived from the cubic param¬ 
agnetic structur d^^ * ^^! The current model of anisotropy change following ultrafast optical 
excitation in this compound assumes a partial reorientation of the {112} spin population, 
either along (111) or {llO}, depending on specific excitation parameterd^. 

In order to disentangle influences arising from the complex domain structure, we prepared 
three different NiO specimens with (111) surface orientation that underwent different anneal¬ 
ing and polishing processes. We prepared samples of 36, 45 and 50 pm thickness. The latter 
two showed bright green color and structural domain sizes of 10-100 pm as well as single 
spin domains, which was verihed using birefringence microscopy and SHG measurementdi^. 
The 36 pm sample showed evenly distributed structural domains smaller than 10 pm and 
was of brownish green color hinting at a deviation from ideal stoichiometrj^i^. 


II. TIME-RESOLVED X-RAY DIFFRACTION 

In order to quantify the response of the magnetic order to optical excitation, we used time- 
resolved non-resonant magnetic x-ray diffraction to measure the sublattice magnetization as 
a function of time after absorption of a femtosecond optical pulse. AFMs like NiO are ideal 
systems to be studied with magnetic diffraction, as the alternating magnetic moments form a 
sublattice with half the periodicity of the structural lattice, giving rise to additional solitary 
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magnetic reflections in reciprocal spacd^. 

The intensity of these magnetic diffraction peaks is proportional to the square of the 
magnetic structure factor (Fm(Q)) where Q is the momentum transfer of the diffraction 
peak!^. 

ihu 

{Fmm = -ro -^(M^(Q)). (1) 

Here, rUe and tq are the mass and classical radius of the electron, c is the speed of light and 
hu gives the photon energy. The magnetic scattering amplitude Mm{Q) is given by 

{M(Q)) = tL(Q).A + S(Q)-B (2) 

where A and B are geometric factors that depend on the scattering geometry and on po¬ 
larization for both the incident and scattered x-rays, S(Q) is the Fourier component at Q 
of the spin density, and L(Q) is the Fourier component at Q of a function related to orbital 
contributions to the magnet izationf^. 

In 1972 Brunei and Bergevin used NiO to demonstrate for the first time the feasibil¬ 
ity of magnetic x-ray diffract iorP^. Since then, the brightness of light sources increased 
tremendously, but the low efficiency of magnetic scattering still impedes time-resolved mea¬ 
surements with ps resolution at conventional light sources. For this reason the experiment 
was carried out at beamline 3 of the x-ray free electron laser (XFEL) SACLA, JaparP^. 

The 45 pm thick (lll)-oriented NiO crystal was mounted on a multi-axis diffractometer 
in horizontal scattering geometry as depicted in Fig. [^a) (tt incident x-ray geometry). To 
measure the intensity of magnetic diffraction from the (| | |) planes, x-ray pulses from the 
XFEL were set to an average photon energy of 7.2 keV, with a full-width-at-half-maximum 
(FWHM) pulse duration below 50 f^^QEII estimated spectral bandwidth of 50 eV. 

At a repetition rate of 10 Hz the total pulse energy was approximately 180 pJ. The Bragg 
angle for the (| | |) planes at this photon energy is 32.4°. The x-ray footprint on the 
crystal was (0.40x0.75) mm^. Polarization analysis using a magnesium oxide analyzer crystal 
in vra-conhguration efficiently suppressed the charge-scattered background, allowing better 
measurement of the the weak magnetic signal. In this configuration the (| | |) reflection 
gives the strongest signal relative to other magnetic reflections, yielding 0.9 photons per 
XFEL-pulse, as shown in the inset of Fig. 

To test for possible contributions from the XFEL second harmonic diffracting from the 
(333) structural planes, we inserted a Si hlter with a nominal thickness of 100 pm into 
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the beam path. This thickness of Si transmits 14% of x-ray radiation at 7.2 keV, bnt 
for any potential second harmonic contribntions at 14.4 keV the transmission is 77%. We 
observed a transmission of the diffracted signal of (12.2 ± 0.6)%, indicating that second 
harmonic contribntions from the (333) strncturally allowed reflection are not signihcant in 
this experiment. 



FIG. 1. Panel (a) depicts the experimental diffraction setup of the horizontal scattering geometry 
at beamline 3 of the XFEL SACLA. An MgO crystal at a Bragg angle of 45° was used to select 
outgoing a radiation that was detected using an yttrium aluminium perovskite detector (YAP). 
Panel (b) shows the spin orientation in the antiferromagnetic phase of NiO for a single spin domain 
of a (lll)-cut crystal where the spin points along |ll2|. 

To excite the sample, an amplihed Ti-sapphire laser system synchronized to the XFEL 
provided pnmp pulses centered at 800 nm wavelength with 50 fs FWHM duration. The 
timing jitter of less than 1 ps exceeded our requirements in order to resolve a potential 20 
ps oscillation period, corresponding to a 50 GHz respons^. The pump pulses were focused 
onto a (0.7x1.7) mm^ (FWHM) spot on the sample at an incidence angle of 24.4°, which 
is close to Brewster’s angle and 8° degrees from the x-ray beam. The incident excitation 
fluence was 37 mJ/cm^, leading to an excitation density of 0.5 ■ 10^° cm“^ near the sample 
surface. As can be seen in section 4, these pump conditions lead to large magnitude SHG 
dynamics, which were interpreted as a reorientation of spins in the [111] directiorP. Spatial 
and temporal overlap of the pump and probe pulses were verihed by measuring the x- 
ray induced optical transmission changes of a GaAs wafer temporarily inserted into the 
sample positiorP^. The intensity attenuation depth was 24 pm and 26 pm for the x-rays 
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and the pump beam, respectively. We measured the intensity of the (| | |) peak at room 
temperature as a function of relative pump-probe delay time over a range of —30 to 30 ps. 
As shown in Fig. 2, the measured relative changes in the diffraction intensity did not exceed 
a value of two times the mean error of photon counting statistics of 2% over the hrst 30 ps 
after excitation. 



FIG. 2. The time resolved diffraction measurements on NiO are shown as blue dots. Broken blue 
lines indicate two times the mean standard error of photon counting statistics, approximately 2%. 
We did not observe a relative change in diffraction intensity exceeding this value. Temporal and 
spatial overlap were verified using an x-ray pump / optical transmission probe experiment on GaAs 
depicted as orange triangles. The inset shows a scan of the NiO Bragg angle 0 in the vicinity of 
the (| I |) refiection. The time scan was carried out at its peak position of 0 = 32.4°. 


III. DISCUSSION OF THE DIFFRACTION RESULTS 

The absence of a change in magnetic diffraction intensity exceeding 2% stands in contrast 
to the SHG response where a drop of more than 20% is seen with similar excitation conditions 
(see Fig. for comparison). The SHG intensity l 2 uj should in principle obey the relationship 
oc where S is the magnitude of the sublattice magnetizatiorP. Based on this one could 
infer that the magnetic order in the excited state drops on the order of several percent. To 
better understand this apparent discrepancy we need to examine more closely what precisely 
the magnetic diffraction signal measures. 

From Eq. [^we see that the magnetic scattering factor is related to a sum of projections 
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of S(Q) and L(Q), Fourier components of the spin density and an orbital density function. 
For our measurement geometry and x-ray polarizations the scattering factor is given by 

Mjro- = 2 sin^ 9 [cos 9 {Li -|- S'!) -|- sin . (3) 


Here, Li and Si are the components of L(Q) and S(Q) along the ( 112 ) direction as dehned 
by Blume and Gibb^and depicted in Fig. [^b). S 3 is the component of S(Q) in the (111) 
direction. We will neglect contributions from Li, although it was found to lead to small 
contributions to the equilibrium magnetic moment in NiCP^. The diffracted intensity is 
then 

I ( 6 ')^o- oc sin^ 9 tan 9 (Si cos 6 * -|- S 3 sin 6 *)^ . (4) 


As discussed in the previous section, the equilibrium sublattice spins can point along any 
of three equivalent { 112 } directions, resulting in three possible spin domains: ( 112 ) (do¬ 
main “A”), (121) (domain “B”), and (211) (domain “C”). Using superscripts to denote 
the different spin domains with their respective ratios of the total population a, satisfying 
= 1 , we have in equilibrium = a^S, Sf’^ = —a^’^S/ 2 , and = 0 . 

The dynamics inferred from previous interpretations of the excited state imply both 
a reduction of the average spin moment and a reorientation of the spin sublattice vector 
along ( 111 ), which in our treatment would lead to a decrease in the Si component and an 
increase in the S 3 component. We can parameterize this change through new time-dependent 
variables a and 7 , with a representing the dimensionless magnitude of the average sublattice 
spin and 7 the reorientation toward the (111) direction. We then have S^ = a a^S cos 
Si’^ = —aO'®’‘"S'cos 7/2 and Ssm.'^. This leads to relative changes in 

the diffraction intensity for each domain type 


m 


( 0 ) 


= (cos 7 -|- tan 9 sin 7 )^ 


m 


( 0 ) 

' B,C 


= a [ - cos 7 — tan 9 sin 7 


( 5 ) 

( 6 ) 


In equilibrium, a = 1 and 7 = 0 . In case of spin reorientation 7 should assume some 
non-zero value. In our data we observe no changes at any time within our estimated 2% 
uncertainty. Since the spin domain population for the probed volume is a priori unknown, 
it is difficult to place a dehnite upper bound for changes in a and 7 . We can, however, point 
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out that the x-ray data show no evidence of strong oscillatory behavior, suggesting that a 
re-examination of the proposed dynamics is warranted. 


IV. TIME-RESOLVED SHG MEASUREMENTS 

SHG has been demonstrated as a powerful tool to study ferroelectricity and antiferro¬ 
magnetism, where the order parameter is correlated to a breaking of the inversion symmetry. 
It is, however, remarkable that the method was proven to be applicable also to centrosym- 
metric materials such as NiO and CoO when the generation is resonantly enhanced and a 
magnetic dipole transition participates in the excitatiorP. In NiO this means, however, that 
SHG is only efficient when using probe light of a narrow band around 1200 nm. Even then, 
the generated second harmonic intensity is very weak, requiring intense probe light pulses 
and long acquisition times. All previous studies were therefore carried out using a funda¬ 
mental wavelength A = 1200 nm and detected the SHG that was emitted in an anti-collinear 
geometry, leading to the observation of oscillations in the SHG response with a frequency of 
either 55 GHz or 1 GH#HZ1. 

Here, we extend these SHG studies with time-resolved measurements in three major re¬ 
spects in an effort to understand these dynamics more completely. First, we carried out 
measurements in both reflection and transmission geometries to test the sensitivity of the 
dynamics to the phase-matching conditions of the SHG process. Second, we performed 
measurements at multiple probe wavelengths to determine whether the observed oscilla¬ 
tions change with probe frequency. Third, we measured the dynamics with different sample 
thicknesses to determine any possible role of the distance between the interfaces. We also 
re-investigated the pump intensity dependence of the dynamics. Our results show that the 
observed dynamics depend strongly on probe wavelength, geometry and acoustic proper¬ 
ties of the samples and appear to be consistent with simulations of the time-dependent 
SHG process when allowing for a strain-induced perturbation of the linear and non-linear 
susceptibilities. 

For these experiments, we employed an amplihed 800 nm Ti-sapphire laser system pro¬ 
viding 100 fs pulses at 1 kHz repetition rate. A fraction of the output beam was fed into an 
optical-parametric-ampliher that generated probe pulses tunable between 1140 and 1400 nm 
wavelength. To allow for the investigation of dynamics extending over several nanoseconds. 


the probe beam made two round trips over a 1.3 m long delay stage. Pump and probe 
beams were combined at a relative angle of 2° using a dichroic mirror. The 800 nm pump 
and IR probe beam were then delivered to the sample with typical pulse energies of 120 
pj and 20 pJ, respectively, and hltered for second harmonic light generated by the optical 
components. After focusing this leads to considerable pump and probe peak fluences of 
50 and 240 mJ/cm^. After the sample, a series of polarizers and color hlters hltered out 
both the fundamental probe and scattered 800 nm light, as well as third harmonic light and 
multi-photon huorescence. The SHG was detected using a GaAs photomultiplier tube. 


A. Reflection geometry measurements 

Using a probe beam wavelength of 1200 nm and detecting the SHG light in reflection 
geometry, we were able to reproduce the SHG response after optical excitation as previously 
reportecPin a 36 pm thick sample that shows evenly distributed small domains, as well as in 
a 50 pm thick sample showing domain sizes of around 100 pm, see Fig. The SHG response 
at 1200 nm features a 55 GHz oscillation, showing an amplitude of about 10% and a decay 
time of above 250 ps. A careful study of the observed oscillation amplitude on fluence and 
excitation pulse length did, however, not show any deviation from linear behavior, as shown 
in the right inset in Fig. This stands in contrast to earlier publications that observed a 
threshold behavior depending on pulse peak powei®®. 

A remarkable observation can be made when detuning the probe wavelength away from 
the resonance condition at 1200 nm. This is experimentally challenging, as the SHG intensity 
is largely reduced, which prevented such detuning in earlier experiments. Nonetheless, over 
an accessible wavelength range between A = 1140 and 1300 nm, we were able to acquire 
data showing a dependence of the observed oscillations on probe wavelength. This suggests 
that the frequencies of the oscillations are not an intrinsic property of the material, but arise 
instead from an interaction between material and the probe beam. These measurements were 
carried out at room temperature and repeated at 150 K to test for temperature dependencies. 
Within the precision of our measurement, we did not, however, observe significant differences 
in the measured frequencies. 

In addition, we observed considerable differences in the SHG yield for different NiO 
samples. SHG from a 1 mm thick bulk crystal was undetectable, while 36 pm and 50 
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FIG. 3. Relative change in SHG in reflection geometry for 3 different wavelengths in the 36 pm 
thick sample. Lines are shown to guide the eye and the measurement of A = 1140 nm is scaled for 
better comparison. Its smaller modulation amplitude is mainly due to an increase in fluorescence 
background for shorter wavelengths. The Fourier transforms and their corresponding peaks are 


shown in the left inset, together with the result for a;+ given by Eq. 16 The right inset shows 
the dependence of the oscillatory amplitude on excitation density at A = 1200 nm (error bars are 
smaller than marker sizes). 


pm thick samples gave observable yield in reflection geometry. As optical properties and 
the domain structure in NiO vary depending on their exact stoichiometrj^^, this might be 
caused by differences in the annealing process. It seems, however, more probable that these 
differences reflect instead a thickness dependence, as the slab thicknesses are of the same 
order as the u and 2u absorption lengths of 23 pm and 35 pm, respectively. This would 
imply that the detected signal in reflection geometry consists largely of internally reflected 
SHG. For this reason, we also studied collinearly generated SHG in a transmission geometry. 


B. Transmission geometry measurements 

Figure 1^ shows the SHG response acquired in transmission geometry over a time inter¬ 
val of 25 ns after excitation for two different slab thicknesses. Using 1200 nm as probe 
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FIG. 4. SHG response to optical excitation in collinear geometry acquired at 1200 nm. The 
oscillations show a recovery in amplitude centered at 10.3 and 13.7 ns for the 36 and 50 pm sample, 
respectively. This is coincident with the observation of phase discontinuities at these times, that 
also occur at 5.15 and 6.85 ns. Pump probe traces show similarity to every other time segment, 
while consecutive segments appear as mirrored (indicated by arrows). 

wavelength, the SHG shows pronounced oscillations with a frequency of about 1 GHz, as 
reported previously for some measurements performed in a reflection geometrjP^. No signs 
of oscillations in the 50 GHz range could be observed. The SHG yield in this conhguration 
is in fact about two orders of magnitude larger than in the reflection geometry. As the 
refractive index for the second harmonic n 2 oo ~ 2.4 implies a Fresnel reflectivity of 17% and 
the absorption length is comparable to the sample thickness for the frequencies considered 
here, this implies that the major part of observed SHG in reflection geometry is in fact due 
to internal reflection of the collinearly created second harmonic light. Moreover, the data 
shown in Fig. [^feature recurrences of the envelope amplitude as well as phase discontinu¬ 
ities that are especially remarkable in the 36 pm thick sample. The observation of these 
phase discontinuities suggests that the recovery of amplitudes is not caused by a beating 
but rather by a reflection. 

There is also a remarkable correlation between the crystal thickness and the observed 
recurrence times. The times at which the phase discontinuities occur corresponds to integer 
multiples of the acoustic round trip time Tac = 2d/vs when using a speed of sound Vs = 
(7.1 ± 0.1) km/s in good agreement with ultrasound measurement^. 

In Fig. I^we also investigated the dependence of the observed dynamics on probe wave¬ 
length in a transmission geometry in the 36 pm thick sample. A pronounced frequency and 
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phase dependence can be observed over a probe wavelength range between 1140 and 1400 
nm. The left inset in Fig. shows frequencies estimated from the Fourier transform based 
on the hrst 5.15 ns after excitation. Phase discontinuities appear for each wavelength at 
5.15 ns, after which the dynamics display mirror-symmetric behavior that manifests as an 
apparent reversal of the direction of the time axis. 



FIG. 5. SHG response to optical excitation in transmission geometry for various probe wavelengths 
in the accessible wavelength range. A common feature of all traces is a time inversion symmetry 
with respect to t = 5.15 ns indicated by a broken line. A strong frequency dependence of the 
induced oscillations is observed. The maximum positions of their Fourier transforms over the first 
5.15 ns is shown in the inset on the left, which also shows the resulting frequencies based on the 
simulation and following Eq. [TB} The fit result of the real part of the refractive index is shown to 
the right. Relevant areas are highlighted. 


V. MODELLING THE SHG RESPONSE TO OPTICAL EXCITATION OF A SPA¬ 
TIALLY CONFINED SEMI-TRANSPARENT CRYSTAL 

For optically absorbing materials, Thomsen et ah presented a model to describe how 
light-induced strain is driven by the thermal stress from a sudden increase in temperature 
after optical excitatiorP^. We will briefly review these hndings in the following. Due to the 
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absorption profile in the crystal and the boundary conditions introduced by the surfaces, 
where applied stress is compensated, the induced stress is strongly inhomogeneous in the 
direction normal to the surface. With density p, strain 7733 , stress (T 33 and the displacement 
in z direction M 3 , (here [HI]), the equations of elasticity can be expressed a^ 

da33 du3 

With the boundary conditions given by the crystal surfaces and the initial condition of zero 
strain before excitation, we obtain solutions 


rimU, t) 


{l-R) 


aQ b 1 + u 
A Cy 1 — V 


F{z,t), 


( 8 ) 


where {1—R) accounts for the reflection loss, with the 800 nm reflectivity R and aQ/A 
gives the excitation density, with a being the 800 nm intensity absorption coefficient, A 
the excited sample area and pulse energy Q. is the material specihc geometric 

response to heat with the expansion coefficient 6 , specific heat capacity Cy and Poisson’s 
ratio u. F{z, t) is the dimensionless solution to the equation of elasticity based on the given 
boundary conditions, given in appendix]^ In our model, an exponential decay term e“*/*'* is 
multiplied with Eq. | 8 ]to approximate the effects of acoustic loss and heat dissipation. Since 
NiO is semi-transparent over the range of relevant wavelengths, it is necessary to consider 
the whole volume of the thin slabs, including strain waves propagating from both interfaces. 
The parameters employed for the calculation are summarized in Table |T} 

As strain implies local variations of interatomic distances, it affects not only the me¬ 
chanical but also the optical properties of the crystal. Typical deformations correspond to 
a relative length change of 10 “^, hence the influence on optical parameters, such as the 
permittivity e, can be treated as perturbations 


Ae{z, t) = 2 (n in) 


dn dn 
+1- 


V33{z,t), 


(9) 


_dr]33 dr]33 

where n and k are the real and imaginary parts of the refractive index. Here we also 
assume that the temperature changes have no direct effect on e. Values of dn/dri 33 7 ^ 0 have 
previously been observed in NiO films as well as in (OOl)-cut bulk NiO using reflectivity and 
ellipticity measurement^^^^. A modulation of the linear refractive index already implies 
possible consequences for the SHG response arising from changes in the SHG coherence 
length. 
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In a similar manner, we can also parameterize possible changes in the second-order sus- 
ceptibity with respect to changes in temperature and strain 




dT 


AT(1 - 

T=RT 


dV33 


V33{z,t). 

»? 33=0 


( 10 ) 


where RT stands for room temperature, denotes the demagnetization time, which is in 
AFM compounds typically few p^, and td is the diffusion time constant td which lies in 
the ns regime. For the sake of simplicity of our model, we do not take into account the 
imaginary component of the nonlinear refractive index. 

The second order susceptibility should depend on the antiferromagnetic order pa¬ 
rameter /, which is related to magnetoelastic lattice distortions^. This in turn depends on 
the temperature, which gives us an expression for the temperature dependence 


x(')(T)cx(1-T/T^)2/^, (11) 

with critical exponent jS = 0.3lP^. The local temperature changes due to optical excitation 
in our experiment are on the order of several K at room temperature, justifying a linear 
expansion in AT. 

A dependence of x^^^ on strain can arise from at least two physical origins. One possibility 
is a strain-induced change in the energies of the {3d)^ states of Ni^"*" that may alter the double¬ 
resonance condition of SHG in NiO. Another possibility is magnetoelastic effects that more 
directly change the sublattice magnetization 1. Both of these possibilities are at present 
beyond our abilities to quantify, and so we simply incorporate 5x*^^V^^33 ^ parameter in 

our model. 

For a full quantitative simulation of the time dependent SHG in NiO, we apply a time- 
resolved propagation matrix calculation. This approach is suitable for our experimental 
geometries as internal reflections are not negligible. To implement temporal and spatial 
variations, the crystal volume is divided along the longitudinal direction in N slices of 
each a few nanometer thickness (well below the scale of optical wavelengths) whose optical 
properties are given by Eq. |^and[^ Restricting the problem to one dimension is permitted 
as the Rayleigh range is many times larger than the crystal thickness. Each slice with index 
m can then be treated as a source of SHG Sm due to the induced nonlinear polarization 
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( 12 ) 


( 2 ) 

where /iq is the vacuum permeability and Xm represents the magnetic dipole assisted non¬ 
linear susceptibility according to Eq. 10 Etirn correspond to the right- and left-ward prop¬ 
agating fundamental helds at slice m. The small value of the second order susceptibility 
allows us to describe the SHG process in the Born approximation in which the fundamental 
helds are independent of the SHG. We can then apply the propagation matrix method 
to derive Et!m at each slice position and calculate the emitted second harmonic helds as 
described in appendix 


Results of the simulation for transmitted and rehected intensities are shown in Fig. and 
respectively. On ns timescales, the SHG in rehection geometry follows the behavior of 
the transmitted SHG, which rehects the fact that its main source is the collinear genera¬ 
tion process and internal rehection, with additional contributions of the SHG arising from 
internally rehected fundamental light. The real part of the refractive index n{u:) used in 
our simulations was obtained from a ht to the experimental data as will be discussed below. 
Further optical parameters employed in these calculations are given in appendixand Table 

El 
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FIG. 6. Simulation of the relative change in transmitted SHG in the presence of optically induced 
propagating strain. The calculation solves the propagation matrix model for a fluence of 50 mJ/cm^ 
and a NiO crystal of 36 pm thickness. Acoustic reflections occur at 5.15 ns and are marked as 
broken lines. The Fourier transform is shown in the inset where the solid line shows the curve given 
by uj- in Eq. 16 The spectrum shows step-like discontinuities due to the acoustic reflections. A 
direct comparison to the data is shown in Fig. 
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FIG. 7. Simulation of the relative change in reflected SHG in the presence of optically induced 
propagating strain. The calculation shown here uses the same parameters as the simulation in 
Fig-i The spectrum is shown in the inset, together with the frequency of Brillouin scattering 


of the fundamental ujb = ‘InVsUjjc and the curve given by a;+ of Eq. 16 These are shown as 


broken and solid lines, respectively. The simulation, in particular for longer wavelengths, contains 
spectral components at ojb- This is not reflected by our measurements and most likely due to 
an underestimation of absorption and overestimation of dn^j/dTj^^ at the fundamental wavelength. 
However, small contributions were previously observed at A = 1200 nnP. 


VI. DISCUSSION OF THE OPTICAL RESULTS 

According to our model, the dynamics in SHG are the result of the time dependent sec¬ 
ond order polarization and light propagation inside the medium in the presence of coherent 
acoustic waves. By using the propagation matrix method the model takes into account 
multiple reflections of fundamental and second harmonic helds as well as the strain waves 
inside the crystal. This is especially important as a major part of the SHG in the reflection 
geometry arises from the internal reflection from the interface on the opposite side of the 
crystal. 

Although to obtain approximate quantitative accuracy the model includes many different 
effects, it is possible to gain some additional insight by considering in general the effect of 
strain waves on phase matching for a collinear or anticollinear SHG process. In a uniform. 
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non-absorbing medium, the output intensity in the hrst Born approximation is related to 
the effective length L by 

2 


hu oc 


X 


e'x.Y){i/S.kz)dz 

( 2 ) 2 

4sin2(AfcL/2) 




(13) 


where Afc = 2 koj^k 2 uj, with the sign depending on whether the fundamental and SHG beams 
are collinear or anti-collinear. Strain-induced modulation of the linear and nonlinear optical 
coefficients causes small, z-dependent changes in both Ak and If we approximate a 

strain wave propagating through the crystal away from the front interface as a real-valued 
discontinuity in both Ak and that moves with the sound velocity Vg, we obtain 


hoj oc 


rvst 


(x'^> + i ( 2 )) exp{i{Ak + 6 Ak)z)dz 


+ 


exp{iAkz)dz 


'Vet 


(14) 


where we assume that beyond the discontinuity the optical constants are unperturbed, and 
before the discontinuity the second-order susceptibility changes by 6 ^( 2 ) and the phase mis¬ 
match Ak changes by Evaluating this to hrst order in 6 ^( 2 ) and d^k yields 


hu oc 



2 p 

Ak 



4sin2(AA;L/2) 


2dAkVs't{s,m{Ak{L — Vst) + sm{Akvst)) 

^v( 2 ) 


X 


( 2 ) 


^ ) ( 2 sm^(AA^L/ 2 ) 


(15) 


-f cos(Afc(L — Vgt)) — cos{Akvst) 

We see immediately from this expression that the intensity of the SHG is modulated in time 
with a frequency 

UJ+/- = Re{ 2 k^ ± k 2 uj)vs. (16) 

where we explicitly take the real part of the expression for Ak to extend our result to the 
more realistic case where there is a small imaginary component to both and k 2 w 


The uj- values from Eq. 16 correspond to the slow oscillations in the transmission ge¬ 
ometry at frequencies near 1 GHz as shown in Fig. and where we show values for 
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U- at various probe frequencies by extracting values of the refractive index from a four- 
parameter Sellmeier equation. The result for n(A) is shown in the right inset in Fig. |^and 
the parameters are given in the appendix in Table This parameterization is only a coarse 
approximation of the linear optical dispersion, since the absorption spectrum of NiO shows 
strong features in the observed spectral range that are not accounted for in this modeP^. 
Both the simulations and Eq. predict that the onset phase of these oscillations should 
depend on the unperturbed value of Ak, which in turn depends sensitively on the probe 
wavelength. Different sample thicknesses L also lead to different onset-phases for the same 
probe wavelength, as observed in Fig. The simulation can also account for the observed 
temporary increase in absolute SHG above the equilibrium level, since the strain modu¬ 
lation is under some circumstances able to effectively improve phase matching for higher 
outcoupling. 


The a;+ values from Eq. 16 correspond to the fast oscillations observed in the reflection 
geometry. These values are plotted in the inset of Fig. and match the measured data quite 

well. 

In order to match the magnitude of the oscillations seen in experiment, the simulations 
require a large magnitude of leading to local changes in of several per¬ 

cent (see Table 0 . The physical origin behind this coupling is unclear, but may be found 
in the same magnetoelastic interaction responsible for the large temperature dependence 
of In principle, our observations could also be brought into agreement with a mag¬ 
netoelastic or flexoelectric contribution to SHG due to d{^) 7 ^ 0 or higher order^. 

The exact dependence on strain, however, is outside the scope of this work as the applied 
experimental techniques do not allow us to distinguish between different possible coupling 
mechanisms. Further insight could be obtained using static SHG imaging of strained crystals 
or by employing a strongly focused probe beam in time resolved SHG in order to disentangle 
contributions arising from different depths inside the crystal. 

The agreement between simulations and experimental data is not exact but they repro¬ 
duce the observed frequencies, lead to effects of similar magnitude and give a quantitative 
explanation for the phase behavior of the oscillations. The relative magnitude of the dy¬ 
namics are in general somewhat underestimated. We found the simulation results to be 
highly sensitive to small changes in the assumed linear optical properties. This sensitivity 
may account for some of these discrepancies. Furthermore, the large probe fluences used in 
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these experiments far exceed the limit of small perturbation and will lead to back-action on 
internally reflected beams which is not taken into account. 

As an alternative explanation for the observed effects, in principle a strain-induced modu¬ 
lation of n also leads to coherent Brillouin scattering of the second harmonic light according 
to ojosc = ‘^k 2 ujVs- For wavelengths closer to the band-gap, the strain dependence of the 
refractive index dn/drj strongly increased, which in this case could lead to a pronounced 
visibility of Brillouin scattering at around 600 nm as compared to the fundamental light 
at 1200 nm. In the simulation, both dn/drj 7 ^ 0 and dx^‘^'*/dr] 7 ^ 0 by themselves can lead 
to modulations at the observed frequencies and also to a small spectral contribution at the 
Brillouin scattering frequency of the fundamental beam. However, using our model, linear 
Brillouin scattering alone can not explain the large magnitude of the observed modulation 
in SHG of up to 30% (see Fig. as it would also lead to similarly strong modulations of 
the fundamental light, which was not observed. 

As a hnal remark it might be surprising that despite the seemingly general nature of a 
strain induced change of the nonlinear susceptibility, the dynamic effects observed in NiO 
were not found in similar compounds such as CoO and KNiFa which also rely on magnetic 
dipole assisted SHG. This may reflect a difference in the magnetoelastic interaction in these 
systems. We note, however, that in NiO these dynamics are strongly dependent on sample 
thickness, absorption and dispersion in Ak. It may be that only a narrow set of experimental 
parameters leads to similar dynamics. As the SHG process in these compounds is highly 
restricted by resonance conditions, the range of these parameters is quite limited. 


VII. CONCLUSION 

Picosecond-time-resolved non-resonant magnetic x-ray diffraction was employed as a tool 
to study sublattice magnetism and suggests that the dynamics in NiO observed with SHG 
may not directly reflect dynamics of the antiferromagnetic order parameter. The low diffrac¬ 
tion efficiency hindered a more precise measurement of the dynamic change of the structure 
factor but the result excludes a large drop in the spin sublattice magnitude as suggested by 
previous models. 

The extended SHG data presented here supplement previous investigations and reveal 
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some aspects that appear inconsistent with previous explanations for the ultrafast dynamics 
of SHG in NiO. The threshold behavior in the presence and frequency of the SHG oscillation, 
which was one of the main arguments for an interpretation in terms of dynamics of the 
order parametei®®, was not verihed by the present experiments. Our SHG studies give 
evidence for an acoustic origin of the dynamics, which becomes particularly apparent in 
the observation of echoes that depend on the acoustic path-length, as well as the probe 
wavelength dependencies for the two observed frequency regimes. 

The choice of crystal dimensions used in the current, as well as in previous SHG studies 
on NiO requires careful interpretation of the observations as the semi-transparency allows 
for multiple reflections. In particular, a direct connection between ultrafast induced changes 
in and I cannot readily be made. Implementing a strain dependence of the linear and 
nonlinear refractive index in a full calculation of the SHG process in a spatially conhned, 
semi-transparent medium gives a self-consistent explanation of the observations in various 
geometries and over a wide range of timescales in which the two observed frequencies can be 
understood as arising from coherent acoustic perturbation that effectively projects the local 
collinear and anti-collinear phase mismatch in SHG into the time domain, with observed 
frequencies uj+- = Re(2/c^ ± k 2 uj)vs. These conclusions might be tested further by a direct 
measurement of the dependence of on strain along the [111] direction. 
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Appendix A: Surface generated strain in a confined medium 

Equation gives the solution for strain induced by impulsive optical excitation. Its 
dimensionless spatial and temporal evolution F{z,t) for thin crystals is then given by 


F{z,t) = f{z) (^1 - ^f{VsT{t))^ 

■ sgn(^ - Vs{t)T{t)) 

+f{d) - z + VsT{t)) 

■ sgn(d - ^ - VsT{t)), 


- ^fi\z-Vsr{t)\) 

- \fi\d- z -VsT{t)\)^ 


(Al) 


with the initial spatial distribution of elastic stress 


fiz') 


^a(2d-z') ^ 


and the effective time r(f) 


(A2) 


{ tmod{d/vs), if tmod{2d/vs) < d/vg 

(A3) 

—tm.od{d/vs) + d/vg, otherwise. 

This definition of r(f) emphasizes the time inversion symmetry with respect to the acoustic 
reflections occurring at multiples of d/vg. 


P [ J?] B [GPa] b Cy td [ns] 

6 . 81 EI 193.8^3 0.26l^ 4.2^3 O.SdI^ 5.0 

TABLE I. Elastic and thermal parameters of NiO used to calculate strain in the simulation. The 
bulk modulus B and Poisson’s ratio v are based on ultrasound measurements on a crystal in (111) 
surface orientatioiP^. The parameter td gives the diffusion time constant used in the simulation. 
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Appendix B: Propagation matrix based calculation of SHG 


The Born approximation allows us to use the propagation matrix approach^ to indepen¬ 
dently solve for the fundamental helds inside the crystal 


C 1 + \ / 

-^uj,M 1 -p, /,\ I 

= Pm-iBJ 


E: 


E: 


(Bl) 

Here, Pm-i is a 2 x 2 matrix connecting the right- and leftward propagating fields at slice 0 
with the helds in slice M. Fields at arbitrary slice positions can be derived by solving for the 
transmitted and rehected helds -F^jv+i using the boundary conditions E^q = Eq 

and ^"^+1 = 0 . 

Given that the lifetime of light inside the crystal is much shorter than X/vg, the propa¬ 
gation matrices Pnit) can be calculated as 

0 

PM{t) = JJ Pm{t), (B2) 

m=M 

where propagation through a single slice is given by 

1 


Pmit) = 


1 — r 



g-^n™(p Q 


0 




(B3) 


with r = (rim+iit) — rimit))/(rim+iit) + Umit)), at normal incidence and Umit) corresponds to 
the time dependent refractive index of slice m, while A represents the chosen slice thickness. 
Surface boundaries are included by setting no{t) = ?7,jv+i(t) = 1- 

The time dependent solutions of E~Q{t) correspond to the familiar results of coherent 
Brillouin scattering of the fundamental beam. By knowing the fundamental held in time 
and space, it is then possible to calculate the emitted second harmonic light by applying the 


propagation matrix approach to each slice as a source of SHG according to Eq. 

0 


E, 

E. 


2uj,N+l 


2uj,N+l 


P2uj,mit) 


.m=N 
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1 

1 E2u, 0 j 
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N 

k=l 


fc+l 


P2ui,m{t) 


,m=N 



(B4) 


Equation B4 can be solved for the emitted SHG helds E^^^ and ^ 2 ^ q using source terms 


according to Eq. 12, as well as the boundary condition E 2 ^ = 0. 
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With the dehnitions 
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n> 
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/c+1 

P2u),m.{t) 

jn=N 



the solutions for SHG in transmission and reflection geometry are given by 


(B6) 


( ^ ^ - S'®’ 'l 

V ^i..o ) \ -is- ) 


(B6) 


The bandwidth can be taken into account by a convolution of the solutions 
with the fundamental light spectra. 


Appendix C: Optical parameters employed for the simulation 

The absorption of NiO in the visible and near-infrared range is strongly affected by 
impurities such as excess oxygeii^. In order to determine an absorption spectrum for the 
simulation that suits the NiO crystals measured in our experiments, we used a polynomial 
£t to a known spectrumP^ ao(a;) and adjusted for the specihc impurity concentration by 
htting a{oj) = aao{oj) b to a. set of three direct absorption measurements that we carried 
out for A = 600, 800 and 1200 nm. These measurements also yielded a real part of the 
refractive index of about n = 2.35 ± 0.05. Due to the limited surface quality of our samples, 
we were not able to measure the dispersion in a static experiment with sufficient precision 
to predict the observed frequencies in the SHG response. 
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Bi 

B2 

Cl [pm^] 

C 2 [pm^] 


1 07733 

1.22 

0.1 

0.01 

0.152 

-1 

-10^ 


TABLE II. Optical parameters used in the simulation. The Sellmeier coefficients are given in Bi to 
C 2 - The strain dependencies are shown normalized to the corresponding equilibrium quantity. The 
value of (^^^)/no matches experimental observation^^ without taking into account the observed 
wavelength dependenc^^^l. The magnitude of chosen to best fit the data and 
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